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Abstract 



"••-, The universality of the continuum limit and the applicability of renor- 

^ malized perturbation theory are tested in the SU(2) lattice gauge theory by 

JL computing two different non-perturbatively defined running couplings over a 

1^ large range of energies. The lattice data (which were generated on the power- 

1^ ful APE computers at Rome II and DESY) are extrapolated to the continuum 

limit by simulating sequences of lattices with decreasing spacings. Our re- 
sults confirm the expected universality at all energies to a precision of a few 
percent. We find, however, that perturbation theory must be used with care 
when matching different renormalized couplings at high energies. 
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1. Introduction 

The continuum limit of lattice QCD is known to exist to all orders of 
perturbation theory [1]. Up to finite renormalizations the limit is independent 
of the details of the lattice formulation and the resulting Feynman amplitudes 
agree with those computed via dimensional regularization. 

At the non-perturbative level the existence and universality of the contin- 
uum limit still needs to be established. One would also be interested to know 
whether renormalized perturbation theory is indeed an asymptotically correct 
expansion of the full amplitudes at high energies. These questions must even- 
tually be decided by rigorous analyses [2], but until this very difficult task is 
mastered, some insight can be gained through numerical simulations. This is 
not trivial either, because one needs rather precise data to be able to observe 
a variation of the chosen physical quantities in the accessible range of lattice 
spacings. A further difficulty is that renormalization and cutoff effects must 
be carefully disentangled from each other. 

A particularly suitable quantity to consider in this context is the running 
coupling a(q), because a large range of momenta q can be probed in this case 
by applying the finite-size scaling technique of ref.[3]. Computations of Q!(g) 
along these lines have previously been performed in the pure SU(2) and SU(3) 
Yang- Mills theories [4-9]. Taking advantage of the massive computer power 
provided by the APE machines at Rome II and DESY, we have now been able 
to extend these calculations in several directions, and it is our aim here to 
discuss the results of this effort. 

The model considered is the pure SU(2) gauge theory, the simplest case 
where the questions mentioned above can be addressed. Compared to our ear- 
lier work on this theory, we have simulated larger lattices (thus moving closer 
to the continuum limit) and generated significantly larger ensembles of field 
configurations. In some cases a Symanzik improved lattice action has been em- 
ployed. As a result one is able to control the continuum limit at a higher level 
of precision and with greater confidence. Additional opportunities to study 
the universality of the continuum limit and the applicability of renormalized 
perturbation theory are obtained by comparing two independent running cou- 
plings, with significantly different low-energy behaviour and different cutoff 
dependence. 

We assume the reader is familiar with the finite-size technology explained 
in detail in refs.[3-9]. Our notations and some information on the numerical 
simulations are collected in sects. 2 and 3. The extrapolation of the numeri- 
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cal data to the continuum limit is discussed in sect. 4. We then address the 
question of how well perturbation theory describes the evolution of the run- 
ning couplings and the relation between them (sect. 5). A discussion of the 
applicability of bare perturbation theory is also included here. After that we 
are well prepared to compute the running coupling ay^(q) in the MS scheme 
of dimensional regularization at large momenta q given in units of a suitable 
low-energy scale (sect. 6). Conclusions are drawn in sect. 7 and a compilation 
of simulation data is included in appendix A. 



2. Running couplings 

In this paper we shall discuss two running couplings, Q!sf(?) and aq^p(q). 
The subscripts SF and TP stand for "Schrodinger functional" and "twisted 
Polyakov loop", respectively, indicating the physical amplitudes from which 
the couplings are extracted. They have already been introduced in our pre- 
vious papers on the subject [4-9] so that here we only briefly recall their 
definitions. For simplicity they are given in the language of the continuum 
theory, but it is straightforward to pass to the lattice formulation (cf. sect. 3). 

2.1 Preliminaries 

As indicated in sect. 1 attention is restricted to the pure SU(2) Yang- Mills 
theory in four dimensions. The gauge potential A^(a;), /^ = 0,...,3, thus 
takes values in the Lie algebra of SU(2), which we choose to be the linear 
space of all traceless anti-hermitean 2x2 matrices. Space-time is assumed 
to be a hyper-cubical box of size L with boundary conditions to be specified 
below. The action is given by 

S[A] = -^ [ d*xtT{F^,F^,}, (2.1) 

^^0 Jo 

where F^j^{x) denotes the field tensor, 

F^,, = d^A, - d,A^ + [A^, A,], (2.2) 

and ^0 the bare gauge coupling. 



The basic idea of the finite-size method of ref.[3] is to study the properties 
of the system as a function of the box size L. In particular, the couplings 
ctSF(?) and Q!tp(?) are chosen to depend on L but on no other external scale. 
They are, therefore, running with L and the momentum q is accordingly given 

by 

q=llL. (2.3) 

There are many different ways to introduce such couplings. Our choices were 
made to meet a number of practical criteria (cf. ref.[6]). 

For any particular running coupling a{q) we define the associated coupling 
g^(L) = 4:TTa(q) and the Callan-Symanzik /3-function 

P(g) = -L^. (2.4) 

In perturbation theory we have 

CO 

^(9) --fy,Kf\ (2.5) 

with 5o = ll/247r^ and hi = 17/967r'* being the usual universal coefficients. All 
other (three loop and higher order) coefficients depend on the chosen scheme. 
Another quantity of interest is the step scaling function a(s,u) which 
describes the evolution of the running coupling under changes of i by a factor 
s [3]. Explicitly, if the coupling is equal to u at scale L, its value at scale sL 
is determined by eq.(2.4) and we may define a(s, u) = g^(sL). Since a(s, u) is 
just an integrated form of the /3-function it is scheme dependent too. 

2.2 Definition of asF 

In this case we impose periodic boundary conditions in the space directions, 

A^{x + Lk) = A^{x), (2.6) 

and require 

{Cfcfx) at xq = 0, 
r t / ^'-'^ 

C^(x) at Xq = L. 

The index k runs from I to 3 and k denotes the unit vector in direction k. For 
the classical boundary fields C and C" we take 

Cfc(x) = T]T3/iL, 



CUx) = (vr - 7?)r3/^i, (2.8) 

where < 77 < vr is an adjustable parameter and t^ the third Pauli matrix. 
The Schrodinger functional, 

Z= /D[A]e-^[^l, (2.9) 



involves an integration over all gauge fields A with fixed boundary values 
C and C" . It may be interpreted as the (euclidean) propagation kernel for 
going from the initial gauge field configuration C at time sq = to the final 
configuration C" at time xq = L. At small couplings it is dominated by the 
field configuration B with least action. For the boundary fields specified above 
we have 

Boix) = 0, 

Bk(x) = [xoC', + (L-xo)Ck]/L, (2.10) 

and the effective action T = —InZ may be expanded in the series 

T = g^'To + T^+g'oT2 + ... (2.II) 

withro=5o'^[^]- 

Power counting and an explicit one-loop calculation suggest that the ef- 
fective action is renormalizable up to a divergent additive constant. So if we 
differentiate with respect to the boundary fields, 

dV 

r' = 1^. (2.1-2) 

the divergent part is removed and a renormalized coupling ^SF may be defined 
through 

Note that L is the only external scale on which ^SF depends. To one-loop 
order its relation to the running coupling in the MS scheme of dimensional 
regularization is (cisf = 9sf/^''^) 

"Ms = "SF + 0.9433 X (asF)^ + •• • (2.14) 
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where both couplings are evaluated at the same momentum q. The relation 
to the bare lattice coupling is now known to two loops [18] (cf. sect. 5). 

2. 3 Definition of citp 

The definition of citp involves correlation functions of the Polyakov loops 

P,{x,,x^,x,)= Tr Ipe/o'-^^^^^^^^j, (2.15) 

with twisted periodic boundary conditions on the gauge potentials, viz. 

A^{x + Li)) = ^,A^{x)^l. (2.16) 

The unitary "twist matrices" fij, satisfy 

f)o = f)3 = 1, f)if)2 = -^2^1, (2.17) 

and the extra phase in the definition of Pi is included to guarantee that this 
observable is a strictly periodic function of the coordinates sq, X2 and ^3. 

Twisted boundary conditions eliminate degenerate toron configurations 
and hence enable standard perturbative computations. The perturbation ex- 
pansion of the two-point correlation function of Pi starts at order g^, while 
the correlation function of P3 has a connected part of order 1. We may thus 
define a running coupling ^^p(P) through 

-2 , .. _ . J Ax2Ax^{P^{Ll2,X2,x^)P^{Q,Q,Qy) 

g^.p(L) - ^ ^^^^^^ ^^^^^^2, x,,X2)Ps{0, 0, 0)*) ' ^'-'"^ 

where the normalization constant k = 14.459 ... is chosen such that gj.p(L) = 
g^ to lowest order of perturbation theory. Linear divergences affecting the 
numerator and the denominator factorize [10] and cancel in the ratio (2.18) 
which is, therefore, renormalizable and well-defined in the continuum limit. 
This has been shown to be correct in an explicit one-loop calculation [7], 
which also yields the relation 

cij^ = ciTP -0.5584X (ciTp)^ + ••• (2.19) 
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between ajp(q) = gj.p(L)/4:TT and the coupling in the MS scheme of dimen- 
sional regularization. 

At large L the choice of boundary conditions becomes inessential and the 
correlations of Polyakov loops in twisted and ordinary directions tend to the 
same value. We thus expect that aq^p(q) converges to k/iir = 1.151 . . . when 
q ^ 0. The infrared behaviours of citp and cisf are thus completely different 
(the latter diverges exponentially [3]). 



3. Lattice formulation and numerical simulation 

3.1. Lattice definition o/cisf and aq^p 

We consider a hypercubic lattice with spacing a and size L*. The lattice gauge 
field is denoted by U(x, i^). In the case of the Schrodinger functional we impose 
periodic boundary conditions in the spatial directions, while at time xq = 
and xq = L we require 

U(x,k)\^^=o = ex]){aCk}, [/(a;, A;)|^„=l = exp{aC[}, (3.1) 

where C'k and C^ are given by eq.(2.8). The lattice action 

^[f^] = ^E^(^)t^il-^(^)} (3.2) 

involves a summation over all oriented plaquettes p on the lattice with U(p) 
being the parallel transporter around p. The weight w(p) is given by 

/ X _ J Cfigo) if p is a time-like plaquette touching the boundary, /„ „x 
t 1 elsewhere. 

Note that in the present context we do not need to specify the weight of the 
spatial plaquettes at sq = and xq = L since their contribution vanishes for 
the boundary values chosen. 

As discussed in ref.[4] the lattice boundaries give rise to 0(a) lattice arti- 
facts which can, in principle, be cancelled by an appropriate choice of Ct(go). 
In perturbation theory we have 

ct{go) = 1 + 4^^5o' + 4'^5o + • • • (3.4) 



Setting Ct(go) = 1 is referred to as tree-level improvement in this paper. For 
1-loop improvement we include [4] 



4^) = -0.0543, 



and for 2-loop improvement [18] 

4^) = -0.0115 

is taken in addition. 

The Schrodinger functional and the effective action T on the lattice are 
defined through 

Z = exp{-r} = /D[C/]e-^[^l, D[C/] = JJ dC/(a;,//). (3.5) 

At small couplings the integral is dominated by (the gauge orbit of) the clas- 
sical field 

V(x,i2) = exp{aBi^(x)} 

with B^{x) given by eq.(2.10). The corresponding action is 

L^ . 2 f (tt - 2?7)a2 



a^ 



ro = 24- sin^ ^ '' . (3.6) 



The running coupling g^Yi-^) i^ ^^"^ defined by eq.(2.13) with 77- derivatives of 
the lattice quantities on the right hand side. 

In numerical simulations the coupling is obtained through 



.2 fT^ _ T./ /^^\ 



aiAD = U ( ^ ) • (3.7) 



Explicitly the "observable" dS/dr] is given by 



-g- = -ciiao)-^ E E i^^ W + ^'^ W} ' (3-8) 



where 



EkU) = -^tr{T3V(x,k)U(x + ak,0)U(x + a6,k)-'^U(x,0)-'^] (3.9) 



(E'l^ is defined analogously at xq = L). 

In the case of the coupling gjp(L) we impose twisted periodic boundary 
conditions, i.e. we initially require U(x + LD^fi) = fl^U(x^ fj,)fll. By a change 
of variables the theory is then rewritten as a theory with ordinary periodic 
boundary conditions, but with a modified Wilson action where the sign of some 
plaquettes is reversed (see e.g. ref.[ll]). In this formulation the Polyakov loops 
are given by 

Pi{xo,X2,xs)= Tr I II C/(a;,l)le-'-^^/^, 
P3{xo,xr,X2)=TT \l[U{x,3)\, (3.10) 

U3 = J 

and gj.p(L) is now again defined through eq.(2.18), where the integrals are to 
be replaced by the corresponding sums over lattice points. 

It should be noted at this point that the constant k is taken to be the 
same as in the continuum theory. On the lattice gjp(L) is, therefore, not 
exactly equal to g^ at tree level of perturbation theory. The coupling is still 
proportional to g^ with a proportionality constant differing from 1 by terms 
of order a^ . 

To reduce statistical fiuctuations in numerical simulations we made use of 
the 1-link integral technique which replaces each link in the Polyakov loop by 
its average value in the external field produced by the neighbouring links [12]. 
If Pk denotes the Polyakov loop Pk averaged in this way, the coupling gj.p(L) 
was thus obtained by computing the expectation values of the observables 

Oi = a^ Yl PiiL/2,X2,X3)Pii0,0,0y, 

03 = a^ Y. ^3(i/2, xi, X2)P3iO, 0, 0)*, (3.11) 



Xi,X2 



and taking their ratio. 



3.2 Simulation algorithm 

As already mentioned, the numerical simulations have been performed on the 
APE computers at DESY and Rome II. In the most commonly used configu- 
ration these computers consist of a mesh of 128 computing nodes. The local 
memory is 1 Mword so that lattices with L/a < 24 can be simulated on subsets 
of 8 nodes, i.e. in this case there are 16 independent lattices on the machine. 
Larger lattices must be distributed over more nodes and the performance of 
the simulation algorithm is then slightly reduced due to the communication 
losses. 

In all cases a hybrid over-relaxation algorithm was employed where Aqr 
exactly microcanonical over-relaxation sweeps are followed by 1 Fabricius- 
Haan [13,14] heatbath sweep. The high-quality random number generator of 
ref.[15] was used. Aqr should be chosen so as to minimize the autocorrelation 
times of the quantities of interest. We have not tried to optimize this param- 
eter, but found, after some testing, that Aqr = L/2a is in general a good 
choice. 

We now need to be a bit a more specific on the order in which the link 
variables are updated, because it turns out (to our surprise) that the autocor- 
relation times depend on these seemingly marginal details (cf. subsect. 3.3). 
Only the lattices which fit on 8 nodes are considered here. In this case there is 
on each node a fraction of size L X (i/2)"^ of the lattice, where the first factor 
refers to the time coordinate. The nodes operate in a SIMD mode and the 
corresponding sublattices are, therefore, updated simultaneously. 

A first possiblity now is to update the link variables U{x,ii) in the order 
where the innermost loop runs over sq, the next over si, then X2, then ^3 and 
the last loop over the direction ji. This is referred to as SF-updating, since it 
is the default order used in simulations of the Schrodinger functional. 

Alternatively one may sweep through the lattice and update all link vari- 
ables U{x,ii) at the current site x before going to the next lattice point. This 
order has been used for the computation of ^^p(i) and is, therefore, called 
TP-updating. 

At L/a = 24 our program with SF-updating achieves a link update time 
of 34 jis (on each node). This number refers to the average update time in 
the mixture of heatbath and over-relaxation sweeps and includes the commu- 
nication overhead. TP-updating is faster by about 25%. For comparision we 
mention that a similiar program, with the same random number generator, 
runs at a speed of about 3.1 jis per link on a CRAY-YMP. 



Table 1. Relative variance V of dS/drj at gly^L) ~ 3.6 



L/a 


V 


L/a 


V 


5 


0.38 


14 


1.6 


6 


0.50 


16 


1.9 


7 


0.64 


20 


2.5 


8 


0.77 


24 


3.1 


10 


1.0 


28 


3.7 


12 


1.3 


32 


4.5 



3.3 Computational efficiency 

An important criterion for the choice of a particular finite volume renormalized 
coupling is the statistical precision that can be achieved in the numerical sim- 
ulations. With this in mind we here discuss the variances and autocorrelation 
times of the two couplings introduced above. 

The relative variance V of any observable O is defined by 

V = {<0^ > - <0>^) / <0>K (3.12) 

In table 1 the variance of dS/dr] is listed for an approximately fixed value 
of ^gp(i). The behaviour of the variance is similiar at other values of the 
coupling. One notices a significant dependence on the lattice spacing roughly 
proportional to {L/aY'^ . 

In the case of the observable Oi [eq.(3.11)], the variance is nearly inde- 
pendent of the lattice spacing in the range 8 < L/a < 16. Its value slightly 
increases from 1.3 at gjp(L) ~ 8.1 to 1.9 at gjp(L) ~ 2.4. The observable 
O3 has a much smaller variance than Oi and will not be considered any fur- 
ther, because the final error on the coupling is dominated by the error on the 
expectation value of Oi. 

Integrated autocorrelation times Tint were estimated by comparing the 
naive errors with our jacknife errors for large bin sizes. We used the Madras- 
Sokal formula [16], which approximates the required four-point autocorrelation 
function by its disconnected part, to calculate the statistical errors on the auto- 
correlation times. The truncation of the range over which the autocorrelation 
function is effectively summed is given by the bin length in our application. 
We found that this kind of error estimate is consistent in the cases where 
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Fig. 1. Autocorrelation times of dS/drj (open symbols) and d (filled 
symbols) for two difi!"erent values of the renormalized coupling in each case. 

independent runs were carried out. 

In fig. 1 we show the integrated autocorrelation times Tint of dS/dr] (with 
SF-updating) and Oi (with TP-updating). We quote Tint in units of sweeps, 
not discriminating between over-relaxation and heatbath sweeps, since this 
unit is roughly proportional to the CPU-time. The figure shows that the 
autocorrelation times are rapidly growing, with no substantial dependence on 
the value of the renormalized coupling. In the present range of L/a the data 
are described by Tint oc (L/a)^ with z ~ 1.0 and z ~ 1.8 for dS/dr] and Oi, 
respectively. It should be added here that the growth of Nqr with L/a was 
chosen to be somewhat slower for TP-updating than for SF-updating (where 
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Nor was always set to L/2a). 

Taken together the computational effort to calculate the renormalized 
couplings ^gp(i>) and gj.p(L) is in both cases scaling roughly proportionally 
to a power of L/a with an exponent equal to 6 or slightly larger than 6. The 
prefactors are such that an order of magnitude less computer time is required 
to obtain ^sf(-^) ^^ ^ given relative precision and for a given lattice size L/a. 

After the computations had been finished we noted that SF-updating, 
when used to compute gj.p(L), leads to a reduction of the autocorrelation 
time of Oi by a factor of 2. Even though the TP-updating program has 
a smaller link update time, it would thus have been profitable to use SF- 
updating in all cases. We have then also discovered that when TP-updating is 
used to simulate the Schrodinger functional, the autocorrelation time of dS/dr] 
increases by a factor 4 already on small lattices. Other updating schemes, such 
as checkerboard ordering, have also been tried, but none of them proved to be 
more efficient than SF-updating. At present we do not know what the cause of 
this unexpected behaviour is. The experience suggests, however, that it may 
be worthwhile to try simple variants of any given simulation algorithm before 
a large scale simulation is started. 



4. Extrapolation to the continuum limit 

We now proceed to discuss our results on the running couplings ^gp(i>) 
and gj.p(L). Some of the data of refs.[5,9] will be included in the analysis. 
The new data are listed in tables 7-9 in appendix A. 

Our principal aim in this section is to show that the extrapolation of a 
number of physical quantities to the continuum limit is well under control and 
that the resulting continuum amplitudes are universal, i.e. that they do not 
depend on the lattice action employed or on any other detail of the lattice 
definitions of the quantities considered. 
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Table 2. Expansion coefficients 60 and 61 for 'Etp(2,u, a/ L) [eq.(4.3)] 



L/a 


^0 


^1 


4 


-0.039219 


-0.014651 


5 


-0.021339 


-0.007222 


6 


-0.016852 


-0.004016 


7 


-0.011382 


-0.002367 


8 


-0.009344 


-0.001564 



4.1 Step scaling function 

On the lattice the step scaling function a(s,u) (subsect. 2.1) is approximated 
by a function S(s, u, a/ L) which describes the evolution o{g^(L) under changes 
of the scale i by a factor s at fixed /3 = 4/^o (^ subscript SF or TP will be 
attached to T,(s,u,a/L) when the corresponding coupling is concerned). The 
defining equation is 

gHsL) = J:(,s,f(L),a/L) (4.1) 



and we expect that 



a(s,u) = lim T,(s,u, a/ L) (4.2) 



(if the continuum limit exists). Simulation results on S(s, u, a/L) are available 
for s = 2 and several values of m = g^(L) and a/L. Our task is to extrapolate 
these data to the continuum limit. 

In perturbation theory the cutoff dependence of the step scaling function 
can be studied by expanding 

T.{2,u,a/L)/a{2,u)= 1 + So{a/ L) + Si{a/ L)u + S2{a/ L)u^ + ... (4.3) 

In the case of Ssf(2, m, a/i) the leading coefficient ^0 vanishes and the 1- 
loop and 2-loop coefficients ^1 and S2 are tabulated in ref.[18]. The expansion 
coefficients of Stp(2, m, a/i) have been computed to 1-loop order (table 2). 

It should be noted at this point that Ssf(2, m, a/i) and the associated 
coefficients Sk have an implicit dependence on the degree of improvement of 
the action. For the tree-level improved action Sk is of order a/L times a 
polynomial in ln(a/L) of degree k — 1. 1-loop improvement implies a reduction 
of the degree by 1 and removes the terms of order a/L from ^1 [which then 
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Fig. 2. Extrapolation of SsF(2,M,a/i) to tlie continuum limit at 
u = 2.037, 2.38, 2.84 and 3.55 (from bottom to top). The leftmost points 
are obtained by linear extrapolation and are slightly set apart for better 
readability. The action is tree-level improved (circles) or 1-loop improved 
(filled circles). 

becomes of order (a/i)^]. And if one employs the 2-loop improved action, the 
order a/L contributions to S2 are cancelled in addition. 

Perturbation theory thus suggests that in the case of Ssf(2, m, a/i) the 
leading cutoff effects are of order a/L. Moreover, the use of an improved 
action should speed up the convergence of Ssf(2, m, a/i) to the continuum 
limit. Our data are in fact compatible with a linear extrapolation in a/L and 
one does observe that the slopes of the lines are reduced when going from 
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Fig. 3. Extrapolation of Ssf(2, m, a/i) (circles) and Sgj, (2, m, a/i) 
(filled circles) to the continuum limit at the same couplings u as in fig. 2. 
The simulations were done using the tree-level improved action. 

the tree-level to the 1-loop improved action (see fig. 2). The plot also shows 
that the continuum limit is independent of the lattice action employed, thus 
confirming the expected universality of the step scaling function (Tsf(2,m). 
Instead of Ssf(2, u, a/L) we may also try to extrapolate 



.(2) 



i:'^{2,u,a/L) 



i;sF(2,M,a/i) 



1 + Si{a/L)u + S2{a/L)u'^ 



(4.4) 



The continuum limits of Ssf(2, m, a/i) and Sgp(2, m, a/i) are trivially the 
same. We may, however, hope that the limit is reached faster in the case of 
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Fig. 4. Lattice spacing dependence of 'Etp(2,u, a/ L) (circles) and 
J:i^^{2,u,a/L) (filled circles) at m = 2.075, 2.430, 3.068, 4.108 and 5.597 
(from bottom to top). The leftmost points are obtained by linear extrap- 
olation. 



Sgp(2, M, a/i), since the cutoff effects are here of order u'^ . Note that this sort 
of improvement can be done with or without improving the action. 

(2) 

The extrapolations of Ssf(2, m, a/L) and Sgp (2, m, a/L) (both calculated 
using the tree-level improved action) are compared in fig. 3. At the largest 
values of the coupling the improvement is evident directly from the figure, 
while the correction seems to go in the wrong direction at the lowest coupling 
u = 2.037. This optical impression is, however, misleading. When the data 
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are fitted by 



s[;p?(2,M,a/i) = asF(2,u) + d(u)a/L, (4.5) 



the coefficients d(u) turn out to be zero wittiin errors for aii vaiues of u. 
In tlie case of Ssf(2, m, a/i), on tlie otiier Iiand, tlie same fit yields slopes 
significantly different from zero for tlie tliree larger values of u. We conclude 
tliat perturbative improvement, i.e. passing from Ssf to Sgp, is, in tliis case, 
very efficient in reducing tlie lattice artifacts. 

In tlie case of Stp(2, u, a/ L) the cutoff effects are of order (a/ L)"^ (cf. ta- 
ble 2). Tlie data for Stp(2, u, a/L) are in fact nearly independent of tlie lattice 
spacing and can easily be extrapolated to tlie continuum limit (see fig. 4). Tlie 
"improved" quantity 

SVY 2, u, a L) = , /7 ' ' ' \^^ 4.6 

^P^ ' ^ 1 + So(a/L) + Si(a/L)u ^ ' 

(wliicli is also plotted in fig. 4) leads to a visible improvement only at tlie low- 
est two values of u. Significant cutoff effects are, Iiowever, observed at larger 
couplings, indicating tliat 2-Ioop or non-perturbative contributions are impor- 
tant for tlie description of lattice artifacts at values of u greater tlian about 
3. Tliis is not really qualitatively different from tlie case of SsF, wliere tlie 2- 
loop terms in eq. (4.4) become important at sucli large values of tlie coupling. 
It is interesting to note in tliis connection tliat tlie cutoff effects in tlie two- 
dimensional non-linear (T-model are also not well described by perturbation 
tlieory unless tlie coupling is very small [3]. 

We conclude tliat perturbative improvement must be applied witli cau- 
tion. No general recommendation can be issued at tliis moment and simula- 
tions of sequences of lattices witli decreasing spacings remain mandatory for 
a safe extrapolation to tlie continuum limit. 

Our final results for tlie step scaling function (Tsf(2,m), togetlier witli a 
description of tlie extrapolation used in eacli case, are compiled in table 3. 
Tlie last two lines in tliis table come from different sets of simulations and 
may be combined to yield (Tsf(2,3.55) = 4.98(4). It sliould be empliasized 
tliat all otiier extrapolations described above yield results compatible witli 
tlie numbers quoted in table 3. 

Tlie values of (Ttp(2, u) used later in tliis paper to compute tlie evolution 
of ^^p(i) Iiave been obtained from Stp(2, m, a/i) by linear extrapolation in 
{a/Ly. Tlie numbers Iiave already been tabulated in ref.[9] and so are not 
reproduced Iiere. 
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Table 3. Extrapolation of the step scaling function Ssf(2, u, 



a/L) 



u 


action 


quantity 


i/ a— range 


o-sf(2,m) 


xVdof 


2.037 


1-loop 


V(2) 


5-12 


2.384(14) 


2.5/4 


2.380 


1-loop 


V(2) 


5-12 


2.85(2) 


4.4/4 


2.840 


1-loop 


y(2) 
^SF 


5-16 


3.59(2) 


4.6/6 


3.550 


tree 


SsF 


7-16 


4.97(6) 


2.8/4 


3.550 


1-loop 


V(2) 
^SF 


7-16 


4.99(6) 


6.9/4 



4.2 Non-perturbative relation between cisf (ind citp 

In the continuum limit there exists a function (l)(u) such that 

-gUL) = <^{-glAL)). (4.7) 

In particular, for small u we infer from eqs.(2.14) and (2.19) that 

(l)(u) = U + 0.1195X u^ + ..., (4.8) 

but (l)(u) is also well-defined beyond perturbation theory. This function thus 
provides another opportunity to study the approach of the lattice data to the 
continuum limit. 

On the lattice the relation between the couplings is of the form 

fML) = ^{fsF(L),a/L) (4.9) 

with some function ^(u,a/L) which converges to (l)(u) in the limit a/L -^ 0. 
The leading cutoff effects are expected to be of order a/L. In table 4 we list 
the two renormalized couplings for a range of L/a, with the bare coupling 
adjusted in such a way that ^sf(-^) — 2.0778. The last column in this table is 
hence equal to ^(u,a/L) at m ~ 2.0778. 

Within errors we do not in this case observe any cutoff effects beyond 
L/a = 6. Again this is striking evidence for the existence of the continuum 
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Table 4. Cutoff dependence of g^p{L) at fixed gly^L) ~ 2.0778 



/3 


L/a 


9sf(L) 


9tp(L) 


3.4057 


4 


2.0778(14) 


3.064(12) 


3.5504 


6 


2.0778(27) 


3.012(15) 


3.6585 


8 


2.0778(35) 


2.954(11) 


3.7464 


10 


2.0778(39) 


2.944(17) 


3.8150 


12 


2.0778(66) 


2.954(13) 


3.8709 


14 


2.0778(52) 


2.965(26) 


3.9200 


16 


2.0778(54) 


2.975(15) 



g'^-p(L) has been computed using the 1-loop improved action 

limit. Extrapolation of the results for L/a > 8, allowing for a linear depen- 
dence on a/L, yields 

ciTP = 0.2374(26) at agp = 0.16535. (4.10) 

We have here just considered one particular (and rather small) value of the 
coupling, but there is now little doubt that the continuum limit of ^(u,a/L) 
is reached similiarly at other values of u, too. 

4-3 Low-energy scale 

So far we have been exclusively concerned with physical quantities defined in 
finite volume. We now introduce an infinite volume low-energy scale r^ and 
compute the value of some reference box size Lq^ in units of tq. We then show 
that the ratio Lq^ /ro converges to a universal value in the continuum limit. 
This result is quite crucial as it allows us to relate the box sizes at which the 
renormalized couplings have been computed to the physical low-energy scales 
in the theory. 

The scale r^ is defined through the implicit equation 

F{ro)rl = 1.65, (4.11) 

where F(r) denotes the force between heavy quarks at distance r in infinite 
volume. The advantages of this particular definition and the numerical com- 
putation of ro have been discussed in ref.[17]. In nature r^ is approximately 

19 



0.9 



0.85 - 



0.8 



0.75 - 



0.7 



0.65 - 



0.6 



0.55 



0.5 



0.45 



0.4 




0.05 



0.1 



0.15 



0.2 



a/L 



0.25 

SF 





Fig. 5. Extrapolation of tlie ratio Lq /ro to tlie continuum limit. The 
data points are from simulations using the tree-level improved (circles), 1- 
loop improved (filled circles) and 2-loop improved (triangles) action. The 
leftmost point represents the result of a linear extrapolation of the data 
generated with the 2-loop improved action. 

equal to 0.5 fm, but it should be emphasized that ro is, in the present context, 
a purely theoretical reference scale with no precise phenomenological meaning. 
The box size Lq^ has already been introduced in refs.[5,6]. It is implicitly 
defined by 

(4.12) 



^^p(io'*) = 4.765. 



To compute the ratio Lq^ /ro one needs to determine Lq^ /a and ro/a at a given 
value of the bare coupling (3. This has previously been done using the tree- 
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level improved and 1-loop improved action [17]. It has already been observed 
there that the use of a 1-loop improved action reduces the lattice artifacts 
in the ratio Lq^ /tq significantly. We have now extended this analysis to the 
2-loop improved action (cf. table 9 in appendix A). 

Fig. 5 shows the lattice spacing dependence of the ratio Lq^ /tq for the 
three different actions. We again expect that the leading cutoff effects are of 
order a/L and thus extrapolate the data linearly. Using the 1-loop improved 
action one obtains Lq^ /ro = 0.573(37) in the continuum limit [17], while with 
the 2-loop improved action the result is 

Ll'^/ro = 0.546(36). (4.13) 

The data from the tree-level improved action are also consistent with this 
value, but are more difficult to extrapolate due to their stronger dependence 
on the lattice spacing. In any case fig. 5 represents further evidence for the 
universality of the continuum limit and underlines the necessity to study the 
cutoff dependence of the physical quantities of interest. 



5. Applicability of perturbation theory 

Perturbation theory is expected to apply at large momenta where the 
renormalized coupling is small. By comparing with our numerical results, we 
can now check to what extent this is true for the quantities considered. 

5. 1 Evolution of the running coupling 

Using the values of the step scaling function (Tsf(2, u) computed in sect. 4, the 
running coupling Q!sf(?) can be calculated at a sequence of momenta q given in 
units of the reference box size Lq^ (see table 5). The details of this calculation 
are as in ref.[5] and so are not discussed here again. It should however be 
emphasized that as a result of our use of the powerful APE computers the 
running coupling is now obtained with a statistical error of less than 1%. 

Since the box sizes in table 5 range from about 0.27 fm to 0.017 fm (if 
we assume r^ = 0.5 fm) it is perhaps not too surprising that the evolution 
of the coupling is rather well described by the 2-loop approximation to the 
Callan-Symanzik /3-function [5]. More precisely if we solve the approximate 
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Table 5. asF at scales given in units of Lq 



SF 



qLr 


asF(g) 


3-loop fit 


1.000 


0.3792 


0.3792 


1.854 


0.2825(16) 


0.2841(13) 


3.597 


0.2260(14) 


0.2262(13) 


7.148 


0.1894(14) 


0.1878(11) 


14.20 


0.1621(13) 


0.1611(9) 



renormalization group equation 

L^ = bof + brf + bfg' (5.1) 

with the initial condition g^(LQ^) = 4.765 and choose 

bf = 0.35(12)/(47^)^ (5.2) 

an excellent fit of the data is obtained (last column in table 5). The effective 
3-loop term included in eq.(5.1) is rather small — its contribution decreases 
from about 7% at the largest coupling in table 5 to 1.4% at the smallest 
coupling. 

The fit may be taken as an economic analytic representation of our data 
(including errors). In the following we shall use it to interpolate ^sf(-^) ^^" 
tween the values listed in table 5. With little risk of running into uncontrolled 
extrapolation errors one may also rely on the fit to compute Q!sf(?) at mo- 
menta q higher than those covered by the data. 

We finally mention that the evolution of the other coupling gj.p(L) is also 
well accounted for by the 2-loop /3-function (cf. ref.[9]). The box sizes here 
range from approximately 0.2 fm to 0.006 fm, while the (relative) statistical 
errors are larger by as much as a factor of 2 compared to the errors quoted in 
table 5. 
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5.2 Relation between cisf (^nd citp 

The relation between our running couplings has already been discussed in 
subsect. 4.2. In particular, we have there been able to compute the value 
of Q!Tp(?) non-perturbatively at a momentum q implicitly determined by 
Q!gp(g) = 0.16535. Considering table 5 this coupling is rather small and thus 
in the range where we should be able to compute Q!tp(?) using perturbation 
theory. 

In the accessible range of couplings the 1-loop term in the series 

ciTP = cisF + 1-5017 X (cisf)^ + •• • (5.3) 

[which one deduces from eqs.(2.14),(2.19)] is however quite large. As a con- 
sequence a straightforward application of the expansion is unsatisfactory. In 
particular, at cisf = 0.16535 the right hand side of eq.(5.3) evaluates to 0.2064, 
thus only accounting for about half of the correction required to reproduce the 
numerical value citp = 0.2374(26). 

From the above we know that at large momenta the evolution of both 
couplings ciTP and cisf is accurately described by the universal 2-loop approx- 
imation to the Callan-Symanzik /3-function. A perturbative relation between 
ctTp(?) and aQ,Y{sq) (note the shifted scale) may, therefore, be expected to 
work well, provided a suitable scale factor s is taken. For such considerations 
the precise choice of s is not relevant, but only the range in which it lies. An 
intuitive choice to try is to set s to the ratio between the A-parameters, viz. 

s = Asf/Atp = 0.27620(2). (5.4) 

The perturbative relation then becomes 

aTp(g) = asF(sg) + 0([Q!sF(sg)]^), (5.5) 

i.e. the scale factor (5.4) has the effect of eliminating the 1-loop term. 

So let us now evaluate eq.(5.5) at the momentum q where cisf = 0.16535. 
Using the effective 3-loop /3-function discussed in subsect. 5.1, one finds 

asvisq) = 0.2289(6). (5.6) 

This is indeed close to the true value of Q!tp(?) given in eq.(4.10) and the 
remaining difference 

[aTp(g) - asF(sg)]/[asF(sg)] =0.7(3) (5.7) 
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could well be accounted for by the higher order terms in eq.(5.5). 

We thus conclude that perturbation theory correctly describes the relation 
between citp and cisf at large momenta g, provided the matching is done at 
appropriatly shifted scales. We have shown that this is the case at a particular 
value of g, but since both couplings evolve with approximately the same /3- 
function, the conclusion is bound to hold for a wider range of momenta. 

5.3 Relation between cisf (ind the bare coupling 
The bare perturbation expansion 

fsF(L) = go + ci(a/L)g^ + c^ia/Ljg^^ + ■■■ (5-8) 

has recently been worked out to 2-loop order [18] and we now would like to 
compare this result with our numerical data for the running coupling. In the 
continuum limit the bare coupling and the lattice spacing a are simultaneously 
taken to zero in such a way that a is an exponentially vanishing function of 
g^. In this context it is hence consistent to drop all terms of order a^, p > 0, 
contributing to the coefficients Ck(a/L), which then reduce to polynomials in 
ln(L/a). 

The resulting series is, however, of limited value, since it can only be 
applied if g^ ln(L/a) is small. A similar situation is encountered in the con- 
tinuum theory, when ^sf(-^) i^ expanded in powers of the coupling ^sf(-^') ^^ 
scales L' much smaller than L. Now if we choose L' to be proportional to a, 
the two series may be combined to obtain an expansion of ^gp(i') in powers 
of the bare coupling with no large logarithms. In terms of ciq = 6*0 /4vr the 
general form of this expansion is 

asF{s/a) = ao + di{s)al + d2{s)al + . . . (5.9) 

where s is a scale factor to be chosen with care. It should be emphasized that 
in this equation cisf (•§/«) is defined in the continuum theory. In particular, 
the momentum s/a must be given in units of some physical reference scale, 
such as Lq^ or ro, for the relation to be meaningful. 

To make practical use of the series (5.9) the scale factor s should be 
chosen such that the low-order coefficients dk(s) are reasonably small (as far 
as possible). Different ways to achieve this have been discussed in refs.[18,19]. 
As already noted above the precise value of s is unimportant, since small 
differences will be compensated by the associated changes of the 1- and 2- 
loop terms in eq.(5.9). As in the case of the relation between cisf and citp 
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Table 6. Estimates of asF(q) at g = 10/a(2.85) 



ClSF 


method 


0.1098 


1-loop in cio 


0.1115 


2-loop in ciQ 


0.1110 


1-loop in ao 


0.1128 


2-loop in ao 



an obvious possibility is to adjust s so that the 1-loop term in the expansion 
vanishes. This leads to the series 

asF(8.83/a) = ciq + 1.287 X al + O(a^). (5.10) 

Note that with this choice of s the 2-loop term is only 1-2% of the leading 
term in the interesting range of bare couplings. 

The unexpectedly large scale factor s appearing in eq.(5.10) has long been 
a source of worry and was interpreted as a sign for ciq to be a "bad" expansion 
parameter [19,20]. It was then suggested to replace ciq by an "improved" bare 
coupling such as 

ao = ao/P, (5-11) 

for example, where P is the average plaquette (in infinite volume) belonging 
to the given value of ^g [20]. The associated expansion of cisf, 

asF(1.17/a) = ao + 0.951 X al + 0(5^), (5.12) 

then comes with a more comfortable scale factor s (which was again fixed by 
requiring the 1-loop term to vanish). 

We now evaluate these formulae at bare coupling (3 = 2.85. The associated 
lattice spacing is denoted by a(2.85). It is fairly small — the corresponding 
momentum l/a(2.85) is about 8 GeV (if we assume tq = 0.5 fm) — and so we 
expect the perturbation expansion to apply. 

From eqs.(5.10) and (5.12) the coupling Q!sf(?) is obtained at two values 
of the momentum q given in units of a(2.85). Due to the contracting nature of 
the evolution, the relative errors at different scales are not directly comparable. 
We hence evolve all results for cisf to the reference energy q = 10/a(2.85). This 
is done using the 3-loop fit of the evolution of cisf discussed in subsect. 5.1. 
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The error quoted in eq.(5.2) is negligible here and the results are listed in 
table 6. 

To compare with the non-perturbative values of cisf (which were deter- 
mined at scales given in units of ig^) we need to compute the conversion factor 
a(2.85)/ Lq^ . This can be done using the tree-level, 1-loop or 2-loop improved 
action. The required information can be found in refs.[5,17] and in table 9. 
As a result one obtains a(2.85)/ig^ = 0.0834(5), 0.0874(5) and 0.0889(5) in 
the three cases. We take the last value and conservatively include 0.0015 as 
a systematic error on the scale ratio. Note however the trivial fact that a 
continuum limit cannot be taken here in contrast to all scale ratios appearing 
within the finite-size method. 

Using the 3-loop fit of our data the numerical value of the SF-coupling at 
q = 10/a(2.85) can now be quoted as 

asF(g) = 0.1135(8), (5.13) 

where the error is obtained by taking the (absolute) sum of the statistical error 
and the systematic error from the scale ratio. 

We conclude that the 2-loop formula (5.12) reproduces the numerical 
value of cisF within 1 error margin. When using the standard bare coupling 
as an expansion parameter, the 2-loop result is off by about 1.8% (2.5 error 
margins). Both series show no pathology if the scale factor s is determined by 
requiring the 1-loop term to vanish. In particular, without independent nu- 
merical control it would seem difficult to us to argue that one of the expansions 
is better than the other. 



6. Computation of a-^ 

We finally discuss the computation of ay^(q) at large momenta q given 
in units of the physical low-energy scale tq (subsect. 4.3). The conversion to 
the MS scheme is of conceptual importance, because all reference to a finite 
volume disappears in this step and one gets a result, which is directly relevant 
to the theory in infinite volume. 

The computation is straightforward. Through the 3-loop fit of our nu- 
merical results, described in subsect. 5.1, the coupling Q!sf(^) is available at 
all momenta q > 1/Lq^ given in units of 1/Lq^ . Using the conversion factor 
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(4.13) we can easily pass to units of I/tq. Then, at large momenta g, the 
1-loop formula 

"ms(^) = "SF(sg) + 0([Q!SF(sg)]^) (6.1) 

is employed, with the scale factor s given by 

s = Asf/A^s = 0.44567. (6.2) 

This choice is motivated by our experience with the perturbation expansions 
discussed in the preceeding section. 

To illustrate the procedure, we choose q = 20/ro which is roughly equal 
to the momentum l/a(2.85) considered above. One then obtains 

a^^iq) = 0.2075(39)(89) (6.3) 

The second error here is equal to Q!sF(sg)"^ and is quoted to remind us of our 
ignorance of the 2-loop term in the perturbative relation (6.1) between a-^ 
and cisF- The associated coefficient is currently being computed and will be 
known soon [18,21] so that this source of error will be better under control. 
The first error in eq.(6.3) is obtained by combining the statistical errors from 
the scale conversion factor Lq^ /ro and the evolution of cisf in quadrature. 

A more precise result on ay^(q) is obtained at larger momenta, because 
the coupling becomes smaller and has a weaker scale dependence. At g = 
200/ro for example (which corresponds to about 80 GeV in physical units), 
one gets 

«M-s(?) = 0.1288(15)(21) (6.4) 

and thus a total error of only 3% (if we assume the perturbative corrections 
to eq.(6.1) are no larger than quoted). 

We would like to emphasize that in these calculations we do not need to 
assume that a-^ evolves according to perturbation theory. The conversion 
from the SF to the MS scheme, using eq.(6.1), takes place in the last step of 
the computation, at the momentum where the MS coupling is desired. The 
evolution of cisf, on the other hand, is numerically controlled. 
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7. Conclusions 

While our results do not prove that the continuum limit of the SU(2) lat- 
tice gauge theory exists in the expected way, they represent further impressive 
evidence for this to be true. All the numerical data are perfectly compatible 
with the assumption, at an unprecedented level of precision. To see this it 
was however important to study the renormalized couplings (and the other 
quantities of interest) on sequences of lattices with decreasing lattice spacings. 
The necessity of such an expensive procedure is made particularly clear by 
fig. 5, where we discuss the relation between physical scales defined in finite 
and infinite volume. 

With our recursive finite-size technique we are able to follow the running 
couplings from low to rather high momenta. As already discussed in our previ- 
ous work [3-9], the evolution of cisf and citp is well described by perturbation 
theory in the whole range of momenta covered. The perturbative relation be- 
tween the couplings must however be applied with care. Our experience is 
that it is important to do the matching at appropriately shifted scales. A 
good choice of scale factor is the corresponding ratio of A-parameters. 

The same comment also applies when considering the expansion of cisf 
in powers of the bare coupling. Here we have observed that the use of an "im- 
proved" bare coupling a^ [eq.(5.11)] leads to slightly better results, although 
there is no sign from the series alone (which is known to 2-loop order) that 
the standard bare coupling ciq would be a bad choice of expansion parameter 
at the appropriate scale. 

We finally note that at the level of precision reached, the perturbative 
formulae employed to convert from one coupling to another must be worked 
out to 2-loop order. In particular, the error on our results for a-^ at large 
momenta is now dominated by the truncation at 1-loop order of the expansion 
of this coupling in powers of cisf • 
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Table 7. Pairs of running couplings at fixed bare coupling (tree-level improved 
action) 



(i 


Lja 


ahi^) 


ahi^^) 


2.7170 


5 


3.5500(59) 


5.364(28) 


2.7903 


6 


3.5500(66) 


5.320(31) 


2.8610 


7 


3.5500(67) 


5.230(31) 


2.9158 


8 


3.5500(75) 


5.151(27) 


3.0091 


10 


3.5500(77) 


5.160(31) 


3.0767 


12 


3.5500(81) 


5.095(36) 


3.1401 


14 


3.5500(70) 


5.089(49) 


3.1874 


16 


3.5500(80) 


5.111(51) 



Appendix A 

We here collect our new data obtained from simulations of the Schrodinger 
functional on the APE computers. We were first interested to improve on 
the accuracy of the old data for the step scaling function Ssf(2, m, a/i) at 
u = 3.55, the largest coupling considered, using the tree-level improved action 
(table 7). 

We have then computed the step scaling function Ssf(2, m, a/i) at the 
same values of u previously considered, but now using the 1-loop improved 
action. These results are collected in table 8. 

Finally the existing data to determine the conversion factor Lq^ /tq [5,17] 
were complemented by computing the /3 values at which ^sf(-^) — 4.765 as 
a function of the lattice size L/a for the case of the 2-loop improved action 
(table 9). 
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Table 8. Pairs of running couplings at fixed bare coupling (1-loop improved 
action) 



(i 


Lja 


ahi^) 


ahi^^) 


3.4448 


4 


2.0370(22 


) 2.366(8) 


3.5207 


5 


2.0370(20 


) 2.375(6) 


3.5880 


6 


2.0370(20 


) 2.386(9) 


3.6463 


7 


2.0370(21 


) 2.383(9) 


3.6984 


8 


2.0370(22 


) 2.401(10) 


3.7840 


10 


2.0370(21 


) 2.384(11) 


3.8506 


12 


2.0370(28 


) 2.383(11) 


3.1775 


4 


2.3800(30 


) 2.869(9) 


3.2576 


5 


2.3800(27 


) 2.886(9) 


3.3201 


6 


2.3800(26 


) 2.890(12) 


3.3852 


7 


2.3800(29 


) 2.860(12) 


3.4342 


8 


2.3800(30 


) 2.888(13) 


3.5205 


10 


2.3800(26 


) 2.888(14) 


3.5934 


12 


2.3800(29 


) 2.870(14) 


3.0151 


5 


2.8400(27 


) 3.589(10) 


3.0797 


6 


2.8400(32 


) 3.607(14) 


3.1343 


7 


2.8400(28 


) 3.645(21) 


3.1900 


8 


2.8400(37 


) 3.622(21) 


3.2759 


10 


2.8400(35 


) 3.602(17) 


3.3423 


12 


2.8400(40 


) 3.598(20) 


3.4035 


14 


2.8400(50 


) 3.592(24) 


3.4556 


16 


2.8400(48 


) 3.622(26) 


2.7700 


5 


3.5500(46 


) 5.060(24) 


2.8355 


6 


3.5500(51 


) 4.967(25) 


2.8946 


7 


3.5500(43 


) 4.964(25) 


2.9442 


8 


3.5500(54 


) 4.994(20) 


3.0257 


10 


3.5500(58 


) 5.038(42) 


3.0974 


12 


3.5500(86 


) 5.080(44) 


3.1579 


14 


3.5500(79 


) 4.966(41) 


3.2048 


16 


3.5500(82 


) 4.924(50) 
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Table 9. Bare couplings vs. lattice size at gly^L) = 4.765 (2-loop improved 
action) 



L/a (3 L/a (3 

5 2.5600(17) 8 2.7279(27) 

6 2.6256(25) 10 2.8061(24) 

7 2.6807(23) 12 2.8756(30) 



References 

[1] T. Reisz, Nucl. Phys. B318 (1989) 417 

[2] J. Magnen, V. Rivasseau and R. Seneor, Phys. Lett. B283 (1992) 90, 

Commun. Math. Phys. 155 (1993) 325 
[3] M. Liischer, P. Weisz and U. Wolff, Nucl. Phys. B359 (1991) 221 
[4] M. Liischer, R. Narayanan, P. Weisz and U. Wolff, Nucl. Phys. B384 

(1992) 168 
[5] M. Liischer, R. Sommer, P. Weisz and U. Wolff, Nucl. Phys. B389 (1993) 

247 
[6] M. Liischer, R. Narayanan, R. Sommer, P. Weisz and U. Wolff, Nucl. 

Phys. B (Proc. Suppl.) 30 (1993) 139 
[7] G. M. de Divitiis, R. Frezzotti, M. Guagnelli and R. Petronzio, Nucl. 

Phys. B422 (1994) 382 
[8] M. Liischer, R. Sommer, P. Weisz and U. Wolff, Nucl. Phys. B413 (1994) 

481 
[9] G. M. de Divitiis, R. Frezzotti, M. Guagnelli and R. Petronzio, Non- 
perturbative determination of the running coupling constant in quenched 
SU(2), Rome preprint (1994), hep-lat/9407028 
[10] A.M. Polyakov, Nucl. Phys. B164 (1980) 171; V.S. Dotsenko and S.N. 

Vergeles, Nucl. Phys. B169 (1989) 527. 
[11] M. Liischer and P. Weisz, Nucl. Phys. B266 (1986) 309. 
[12] G. Parisi, R. Petronzio and F. Rapuano, Phys. Lett. B128 (1983) 418. 
[13] K. Fabricius and 0. Haan, Phys. Lett. B143 (1984) 459. 
[14] B. J. Pendleton and A. D. Kennedy, Phys. Lett. B156 (1985) 393 

31 



[15] M. Liischer, Comp. Phys. Comm. 79 (1994) 100 

[16] N. Madras and A. D. Sokal, J. Stat. Phys. 50 (1988) 109 

[17] R. Sommer, Nucl. Phys. B411 (1994) 839 

[18] R. Narayanan and U. Wolff, Two-loop computation of a running coupling 

in lattice Yang-Mills theory, to appear 
[19] G. P. Lepage and P. B. Mackenzie, Phys. Rev. D48 (1993) 2250 
[20] G. Parisi, in: High-Energy Physics — 1980, XX. Int. Conf., Madison 

(1980), ed. L. Durand and L. G. Pondrom (American Institute of Physics, 

New York, 1981) 
[21] M. Liischer, A. van de Ven and P. Weisz, work in progress 



32 



